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We study the effects of CP violation on the nature of the chiral transition within the linear sigma 
model with two flavors of quarks. The finite-temperature effective potential containing contributions 
from nontrivial values for the parameter 6 is computed to one loop order and their minima structure 
is analyzed. Motivated by the possibility of observing the formation of CP-odd domains in high- 
energy heavy ion collisions, we also investigate the behavior of the effective potential in the presence 
of a strong magnetic background. We find that the nature of the chiral transition is influenced by 
both 6 and the magnetic field. 



I. INTRODUCTION 

In Quantum Chromodynamics (QCD), the existence 
of instanton configurations for the gauge fields [T] and 
their intimate connection with the axial anomaly [5] al- 
lows for a nontrivial topological term in the action. This 
term is usually neglected, since it breaks the original CP 
symmetry and so far experiments indicate that its coef- 
ficient, known as the 9 parameter, is vanishingly small, 
9 < 10~ 10 [3] (see also [4]). The reason why 9 is so small 
(or zero) is unclear, and this issue is called the strong CP 
problem . 

Although it has been proved that CP can not be spon- 
taneously broken in the vacuum of QCD for 9 = [BJ, this 
theorem might not hold for QCD matter at finite tem- 
perature or density [7j. With this premise, it has been 
proposed in Ref. [8] that hot matter produced in ultra- 
rclativistic heavy ion collisions could exhibit domains of 
metastable states that violate CP. These states could be 
described by a QCD action that incorporates the topo- 
logical 0-term, and would decay via CP-odd processes. 
Possible experimental signatures for the presence of CP- 
odd domains are based on charge separation of hadronic 
matter produced in heavy ion collisions [9]. This effect 
is enhanced by the presence of a strong magnetic back- 
ground in the case of noncentral collisions, as was realized 
more recently, by a mechanism that has been called chiral 
magnetic effect [10 , and could in principle be observed a 
RHIC and the LHC. 

In this paper we investigate how the chiral transition 
is affected when the ingredients mentioned above, i.e. 
CP-odd effects and a strong magnetic background, are 
present, as in the case of noncentral heavy-ion collisions. 
For this purpose, we adopt the linear sigma model cou- 
pled with two flavors of quarks as our effective theory |llj 
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1 Frameworks that preserve the CP symmetry were also proposed, 
the most studied being axion theories [5], but those were not 
confirmed by experiments to date either, leaving the problem 
unsolved. 



to study the chiral transition 2 . To include the effects of 
the presence of the axial anomaly and CP violation, we 
add a term that mimics the presence of nontrivial gauge 
field configurations, the 't Hooft determinant [26 . The 
latter is a function of the parameter 9, and is respon- 
sible for CP violation for non- vanishing values of 9. In 
this framework, the presence of CP violation is directly 
related to a nonzero r/ condensate. 

The rich vacuum structure brought about by a nonzero 
9 term in the action [37] has an influence on the chiral 
transition, and generates a more complex picture in the 
analysis of the phase diagram of strong interactions 3 . As 
will be shown below, three phases emerge: one with a and 
•q condensates, another where the rj condensate vanishes 
and the a condensate remains, and a phase where both 
condensates are melted. The topography of extrema is 
therefore rich, and metastable minima do appear in cer- 
tain situations. 

Once we have an effective theory describing chiral sym- 
metry restoration at finite temperature in the presence 
of CP violation, we study how this CP-odd linear sigma 
model is affected by the presence of a strong magnetic 
background, as is presumably generated in the case of 
noncentral high-energy heavy ion collisions [10 . This 
extends the analysis we presented in Ref. [23], where we 
investigated the effects of a strong magnetic background 
on the nature and dynamics of the chiral phase transi- 
tion at finite temperature and vanishing chemical poten- 
tial, and found that the nature of the chiral transition is 
modified. Here we find that the main result obtained in 
[2"9"] , namely that the presence of a strong magnetic back- 
ground turns the crossover into a weak first-order phase 



2 This effective theory, especially the <r-7r-quark sector, has been 
widely used to describe different aspects of the chiral transition, 
such as thermodynamic properties rT3l fl4l fT5l fl6l fTTl rT51 [Tol 
1201 1211 22 and the nonequilibrium phase conversion process 23 j5 
as well as combined to other models in order to include effects 
from confinement 1241 1251 , usually without the inclusion of the 
term. 

3 See, for instance, Ref. 1281 for a detailed study of the phase 
structure of the two-flavor NJL model at 9 = 7r. 
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transition, is maintained. In fact, the combination of 
effects from high magnetic fields and the CP-odd contri- 
bution allows for nontrivial combinations of metastable 
CP-odd and chirally symmetric phases, a structure that 
might be relevant for the supercooling dynamics that pre- 
sumably happens after a high-energy heavy ion collision. 

Modifications in the vacuum of CP-symmetric QCD 
by the presence of a magnetic field have been investi- 
gated previously within different frameworks, mainly us- 
ing effective models [30l EU E21 (33l EH EU ESI [37], espe- 
cially the NJL model [33], and chiral perturbation theory 
[3UI HOI HI], but also resorting to the quark model [4"2"] 
and certain limits of QCD [53]. Most treatments have 
been concerned with vacuum modifications by the mag- 
netic field, though medium effects were considered in a 
few cases, as e.g. in the study of the stability of quark 
droplets under the influence of a magnetic field at finite 
density and zero temperature, with nontrivial effects on 
the order of the chiral transition [H]. More recently, 
magnetic effects on the dynamical quark mass [45 and on 
the thermal quark-hadron transition |46] , as well as mag- 
netized chiral condensates in a holographic description of 
chiral symmetry breaking |47j . were also considered. 

The paper is organized as follows. Section II presents 
the low-energy effective model adopted in this paper to 
investigate the chiral transition in the presence of CP vi- 
olating terms. In Section III we show our results for the 
effective potential at finite temperature, and discuss the 
condensates. In section IV we incorporate effects com- 
ing from the presence of a strong magnetic background. 
Section V contains our conclusions and outlook. 



II. CP-ODD LINEAR SIGMA MODEL 

The classical Lagrangian for QCD is invariant under 
parity, P, and charge conjugation, C, transformations. 
However, at a quantum level these symmetries are not 
preserved. This breaking manifests itself as the axial (or 
ABJ) anomaly [2J, and the current associated with this 
symmetry is no longer conserved, even in the massless 
limit for quarks: 

where g is the gauge coupling, Nf is the number of quark 
flavors with masses irif , F is the dual of the gauge field 
strength tensor F, and sums over color indices a and 
flavor indices / are implicit. Naively, the last term in 
the expression above would be irrelevant, since it can 
be written as a total derivative. However, if one con- 
siders topologically nontrivial gauge field configurations, 
this term is nonvanishing and can be associated with the 
winding number of each configuration. Each sector of 
configurations, with a given winding number, has its own 
vacuum, and the true vacuum of the theory becomes a su- 
perposition of the several inequivalent topological vacua 
\9) = X)n e ~ m9 l n )> where is a free parameter. In this 



context, to calculate expectation values one has to com- 
pute an average over a sector with fixed winding number 
and then sum over all sectors. This procedure is equiv- 
alent to adding a topological term in the Lagrangian, as 
follows: 

Cqcd = Cd - ^g 2 F^F^ . (2) 

To describe the chiral phase structure of strong inter- 
actions including CP-odd effects, we adopt an effective 
model that reproduces the symmetries of QCD at low en- 
ergy scales, and has the appropriate degrees of freedom at 
each scale: the CP-odd linear sigma model coupled with 
two flavors of quarks. The chiral mesonic sector is built 
including all Lorentz invariant terms allowed by symme- 
try and renormalizability. Following Refs. [26, 48J, one 
can write 

-^Tr[(^)2] + I [ e «> det(0) + e~' w det(0+)] 

+Tr[h(ct> + ^)). (3) 

The potential in the Lagrangian above displays both 
spontaneous and explicit symmetry breaking, the latter 
being implemented by the term ~ h. The strength of 
CP violation is contained in the 't Hooft determinant 
term, which encodes the Levi-Civita structure of the ax- 
ial anomaly and depends on the value of the parameter 
9. 

Expressing the chiral field <j> as 

4>= ^=(a + ir]) + ~(a+in)-T , (4) 

the potential takes the following form (substituting the 
parameter h by H = \/2h): 

V x = ~(a 2 + 7T 2 + V 2 + S§) 

-|cos0 (a 2 + 7? 2 + T] 2 + al) 
+c sin 9 (err/ — 7? • do) — Ha 
+ l(X 1 + ^)(a 2 + V 2 + n 2 + d 2 ) 2 
2A2 

+ —^- (a d + r]ir + 7T x a ) 2 ■ (5) 

The parameters a, c, H , Ai and A 2 can be fixed by vac- 
uum properties of the mesons, so that the effective model 
reproduces correctly the phenomenology of QCD at low 
energies and in the vacuum, for vanishing 9, such as the 
spontaneous (and small explicit) breaking of chiral sym- 
metry and experimentally measured meson masses. 

In our treatment, the chiral mesons are coupled in the 
standard Yukawa fashion to two flavors of quarks. The 
latter constitute a thermalized fluid that provides a back- 
ground in which the long wavelength modes of the chiral 
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fields evolve. The full Lagrangian has the form: 



C = 



1 



1 



1 



- V(a,T}, ff, (S ) + ^f{il^d^ } 

— gtjjf(—ij 5 T-ir + <T'-f 5 T-do — i , ~f 5 r])'tpf 



(6) 



For T — quarks degrees of freedom are absent (excited 
only for T > 0), and results from chiral perturbation 
theory for the broken phase vacuum are reproduced |49j . 
For T > 0, quarks are relevant (fast) degrees of free- 
dom and chiral symmetry is approximately restored in 
the plasma for high enough T. In this case, we incorpo- 
rate quark thermal fluctuations in the effective potential 
for the mesonic sector, i.e. we integrate over quarks to 
one loop. 

To one loop, and within a classical approximation for 
the chiral field background, a standard functional inte- 
gration over the fermions (49] gives an effective potential 
of the form V e ff = V(<fi) + V q (</>), where the quark con- 
tribution is given by 



V„ 



-v q T 



d 3 k 
(2tt) 3 



In 1 + e 



-E k (<P)/T 



(7) 



where v q — 24 is the color-spin-isospin-baryon charge 
degeneracy factor, and E k {$) = (k 2 + M 2 ^)) 1 / 2 , M q 
playing the role of an effective mass for the quarks. 
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Following a mean field analysis, we take a = (a) + a' 
and rj = (r/) +77', and assume that the remaining conden- 
sates vanish. In this approximation, the effective mass for 
the quarks is simply given by M q = g^J (a) 2 + {rj) 2 . The 
effective potential can then be split into a classical piece 
and its fluctuation corrections, besides the contribution 



coming from the quarks, i.e., V e // — , , 
The classical contribution has the form 



yfluct 



v, 



^((a) 2 -v 2 ) 2 -H(a) 



((V) 1 



Ug) 2 — c sin9 (a) (77) 
A 



-(a) 2 ( v ) 2 ~-(vl + 4) 



(8) 



where we have defined the combination of couplings A = 
Ai + A 2 /2, and the following auxiliary quantities: 



,2 



1 cos 



A 



2 2 

Uf, = Va 



2c 



■ cos( 



(9) 



The fluctuation correction, up to quadratic terms, is 
straightforwardly worked out and provides the means to 
fix the parameters of the model to meson masses and 
the pion decay constant in the vacuum in the ab- 
sence of CP violation. Concretely, we impose that in 
the vacuum chiral symmetry is spontaneously broken and 
the expectation value of the chiral condensate is given 
by (er) = /tt = 93 MeV. The coefficient of the term 
that breaks explicitly chiral symmetry is given by the 
PC AC relation H = /^m 2 . All the other parameters are 
fixed to reproduce the mesons masses: m a = 600 MeV, 
m, = 574 MeV, m,, = 138 MeV, and m ao = 980 MeV. 

This completes the setup of the CP-odd linear sigma 
model to be used as an effective theory in the investiga- 
tion of the chiral transition in the presence of nontrivial 
CP violating processes. In what follows we present our 
results for the effective potential for representative values 
of the parameter and the temperature. Effects from a 
strong magnetic background will be incorporated later. 




(c) (d) 

FIG. 1: Contour plots of the effective potential in the vacuum, 
(a) 9 = 0; (b) 6 = tt/4; (c) 6 = tt/2; (d) 9 = w. The numerical 
values are in MeV. 



III. RESULTS FOR THE EFFECTIVE 
POTENTIAL 

The effective potential we obtained in the previous sec- 
tion is a function of the condensates of a and 77, and of 
the CP violation parameter 6. Below we analyze the de- 
pendence of the effective potential on 9 for different val- 
ues of the temperature. In the vacuum, i.e. in the zero 
temperature case, the potential has two minima in the 
a direction for 9 = 0, as expected from usual (CP-even) 
linear sigma model. The contour plot for the effective 
potential in the plane (a) — (77) for this situation is pre- 
sented in Fig. [lja). As we increase 9, the pair of minima 
rotates, and for 9 = n it is almost in the {rj) direction, 
as can be seen in Figs. |T|(b) and (c). Although this is 
not a realistic (physical) case, since nonzero values of 9 
are not observed in the vacuum of strong interactions, it 
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is nevertheless interesting to study all the phase space to 
have a complete map of the dependence on 9. 
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(a) (b) 



FIG. 2: Contour plots of the effective potential for 9 = 0. (a) 
T = 120 MeV; (b) T = 160 MeV. The numerical values are 
in MeV. 



temperature, the minima move towards the center, in- 
dicating a chiral symmetry restoration, as displayed in 
Fig. [4] However, in this case, there is a barrier between 
the global minimum and the new minimum that will be- 
come the true global minimum at high temperature at 
r\ = 0. In contrast to the previous case, this signals a 
first-order transition, and the possibility of metastable 
CP-odd states. 




(a) (b) 



Keeping 9 — and increasing the temperature, the 
two minima remain in the (cr) direction, as expected, and 
both minima approach to the origin, (cr) ~ 0, so that for 
high enough temperatures chiral symmetry is approxi- 
mately restored (see Fig. [2]). Chiral symmetry is not 
completely restored due to the explicit breaking caused 
by the mass term (~ H) in the effective potential. In Fig. 
[3] we plot the potential in the (cr) direction for 9 = and 
different values of the temperature. As expected from the 
usual linear sigma model, the transition to the phase with 
approximately restored chiral symmetry is a crossover. 



T = 


120 MeV 


T = 


140 MeV 


T = 


160 MeV 




180 MeV 



> 




75 100 
<o> 



FIG. 3: Effective potential normalized by the temperature in 
the (a) (in MeV) direction at 6 = for several values of the 
temperature. 



For 9 = 7r, the degenerate minima are almost in the 
(77) direction, as it is shown in Fig. |ljd) . Increasing the 



FIG. 4: Contour plots of the effective potential for = n. (a) 
T = 125 MeV; (b) T = 128 MeV. The numerical values are 
in MeV. 
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FIG. 5: Effective potential normalized by the temperature in 
the (77) (in MeV) direction at = n for different values of the 
temperature. 



To show more clearly the barrier, we display in Fig. [5] 
the effective potential at 9 — tt for different values of the 
temperature. One can notice the new minimum emerging 
at (77) = for a temperature around T — 126 MeV. 
In this case, the critical temperature is lower than the 
critical temperature for the melting of the condensate 
(cr). Since the critical temperatures for the melting of (cr) 
and (77) are different, three different phases are allowed 
in systems with 9 between zero and tt: one in which 
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both condensates are present, another one where the (rj) 
condensate vanishes, and a phase where both condensates 
vanish. 
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FIG. 6: Absolute value of the condensates (in MeV) as func- 
tions of the temperature. As 9 approaches n, the transition 
in the {rj) direction becomes stronger. Here full lines denote 
(a) and dotted lines (rj). 
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FIG. 8: Contour plots of the effective potential for 8 = n 
and different values of B. The temperature was chosen to be 
T = 195 MeV. (a) B = 5ml; (b) B = Wml; (c) B = 15ml; 
(d) B = 20ml. The numerical values are in MeV. 



phase transition. It is also clear in the figure that the 
critical temperatures for the melting of the two conden- 
sates are different, allowing for the possibility of three 
different phases, as mentioned above. 
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IV. EFFECTS FROM A STRONG MAGNETIC 
BACKGROUND 



FIG. 7: Contour plots of the effective potential for 8 = 
and different values of B. The temperature of each plot was 
chosen to be close to the respective critical temperature, (a) 
B — Wml; (b) B = 20mJ. The numerical values are in MeV. 

In Fig. [6j we plot the absolute value of the expecta- 
tion values of a and rj as functions of the temperature 
for different values of 9. Here, full lines stand for the a 
expectation value and dashed lines represent the conden- 
sate (?/) . Since their values are "complementary" , the 
curves for different values of 9 can be paired using the 
value of each condensate at T = 100 MeV: for 9 = 
the highest solid line corresponds to (a) and the lowest 
dashed line to (ij), for 9 = 7r/4 the second to highest solid 
line corresponds to (a) and the second to lowest dashed 
line to (?y), and so on. One can see that while (a) decays 
smoothly for different values of 9, the condensate (rj) de- 
cays sharply when 9 approaches 7r, indicating first-order 



Motivated by the possibility of observing the formation 
of CP-odd domains in high-energy heavy ion collisions in 
the presence of strong magnetic fields for noncentral col- 
lisions, we can incorporate the effect of a strong magnetic 
background on the effective potential for the CP-odd lin- 
ear sigma model, and study how the chiral transition and 
the condensates are modified. Assuming that the system 
is now in the presence of a strong magnetic field back- 
ground that is constant and homogeneous, one can com- 
pute the modified effective potential following the proce- 
dure we presented in detail in Ref. |29j . In what follows, 
we simply sketch the method and move to the discussion 
of the results. 

For definiteness, let us take the direction of the mag- 
netic field as the z-direction, B = Bz. The effective 
potential can be generalized to this case by a simple 
redefinition of the dispersion relations of the fields in 
the presence of B 7 using the minimal coupling shift in 
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latter provide the new dispersion relations 
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FIG. 9: Effective potential normalized by the temperature 
for 6 — n in the (rj) (in MeV) direction. The plots, from 
the top to the bottom of the figure, correspond to B — 5ml, 
B = 10ml, and B = 20ml. 



the gradient and the field equations of motion. For 
this purpose, it is convenient to choose the gauge such 
that A» = (A , A) = (0, -By, 0,0). Decomposing the 
fields into their Fourier modes, one arrives at eigenvalue 
equations which have the same form as the Schrodingcr 
equation for a harmonic oscillator potential, whose eigen- 
modes correspond to the well-known Landau levels. The 



Pin = pl+m 2 + (2n+l)\q\B: 



Pan 



p~ 



(2n + 1 -<j)\q\B. 



(10) 
(11) 



for scalars and fermions, respectively, n being an integer, 
q the electric charge, and a the sign of the spin, ft is also 
straightforward to show that integrals over four momenta 
and thermal sum-integrals acquire the following forms, 
respectively: 



T E 



d 4 fc 

(2tt)4 

d 3 k 
(2^)3 



2tt 



n=0 
\q\BT 



£/'f^, (12) 

* — ' ' Z7T Z7T 



2tt 



EE 

l n=0 



dk z 
2ir 



(13) 



where n represents the different Landau levels and I 
stands for the Matsubara frequency indices [IS] ■ 

In our effective model, the vacuum piece of the poten- 
tial will be modified by the magnetic field through the 
coupling of the field to charged pions. To one loop, and 
in the limit of high B, eB 3> ml, one obtains (ignoring 
contributions independent of the condensates) [22] 



vY + 



VY = 



2mieB 



32tt 



2 



log 2. 



(14) 



Thermal corrections are provided by pions and quarks. 
However, the pion thermal contribution as well as part 
of the quark thermal contribution are exponentially sup- 
pressed for high magnetic fields, as has been shown in 
Ref. [55] . The only part of the quark thermal part that 
contributes is 



V q T = -N t 



eBT 2 
2tt 2 



dx In f 1 + e 



-y/x 2 +M%/T* 



(15) 

where N c = 3 is the number of colors. Therefore, the 
effective potential computed in the previous section is 
corrected by the contributions in (14) and (15) in the 



presence of a strong homogeneous magnetic background. 

For 9 = 0, i.e. in the case of the CP-even linear 
sigma model, the situation is the same as the one we 
investigated in Ref. [29]. Initially, there are two min- 
ima in the (a) direction, one global and one local. As 
the temperature increases, they move towards the cen- 
ter, approximately restoring chiral symmetry (not fully 
restored since there is an explicit breaking term in the ef- 
fective potential, due to the nonzero quark masses). For 
large enough values of the magnetic background, the na- 
ture of the phase transition is altered, becoming of first 
order instead of a crossover, in accordance with the find- 
ings of Ref. [59] • In Fig- [7] we show contour plots of 
the full effective potential at 9 — for different values 
of the magnetic field. In each plot we display the poten- 
tial close to the critical temperature, so that the case of 
B = 10ml is plotted at T = 120 MeV, while the case of 
B = 20ml is plotted at T = 100 MeV. For B = 5m£, the 
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barrier is quite small and the critical temperature higher, 
as discussed in Ref. [25] . 

Augmenting the value of to 7r, we can study the be- 
havior of the rj condensate. As in the case in the ab- 
sence of a magnetic background, this condensate under- 
goes a first-order phase transition. Moreover, in the pres- 
ence of a strong magnetic field, its critical temperature 
is strongly affected. Analogously to our findings in Ref. 
[2"9"] for the CP-even case, a field B = 5m% increases the 
critical temperature, but now for the rj condensate. The 
behavior of the critical temperature for stronger values 
of the field is analogous to the one of (a), i.e. it drops 
considerably In Fig. [8] we plot the effective potential at 
temperature T = 195 MeV. For B = 5mJ, this temper- 
ature is not enough to take the system to the critical re- 
gion. On the other hand, for B = IOto^ this temperature 
is close to the critical temperature, and for B = 20m^ it 
is larger than the critical value. This can also be clearly 
seen in Fig. [9j where we display slices of the effective po- 
tential in the rj direction for different values of magnetic 
field and temperature. 

Besides modifying the value of T c , the presence of a 
strong magnetic background deepens the absolute mini- 
mum and favors a first-order transition. 



V. CONCLUSIONS AND OUTLOOK 

The possibility of probing the vacuum structure of 
QCD in high-energy heavy ion collisions, bringing some 
light to the understanding of the strong CP problem, is 
a very exciting prospect. The consideration of the theo- 
retical description of the environment that might be pro- 
duced under the appropriate experimental conditions, es- 
pecially the major role played by the presence of a strong 
magnetic background, has already brought up very inter- 
esting phenomena, such as the chiral magnetic effect [TU] 
and the possibility of converting the nature of the chiral 
transition from a crossover to a first-order phase transi- 
tion gpj. 

The complete physical scenario in heavy ion collisions 
is rather complicated [9j [10], so that it is more pru- 
dent to consider, theoretically, the role of each relevant 
ingredient separately, aiming at a more complete picture 
at the end. Previously, we have investigated the effect of 
a strong magnetic background on the nature of the chi- 
ral transition |29j . obtaining remarkable effects as men- 
tioned above, and opening a new line of possibilities in 
the study of the phase structure of strong interactions. 
In this paper, we investigated how the chiral transition, 
more specifically how the effective potential written in 
terms of the condensates (a) and (17), is altered when 



one includes the possibility of strong CP violation and a 
strong magnetic background. For this purpose we built 
a CP-odd extension of the linear sigma model coupled 
with two flavors of quarks at finite temperature and zero 
density. The inclusion of magnetic effects is implemented 
via a redefinition of the dispersion relations. 

We found that, in the absence of magnetic fields, 
the a condensate behaves pretty much like it does in 
the usual CP-even linear sigma model, changing from 
a phase where chiral symmetry is broken ((cr) 7^ 0) to a 
phase where it is approximately restored ((cr) ss 0) via 
a crossover as the temperature is increased. For non- 
vanishing values of the parameter 6, a condensate of the 
field 77 builds up for low temperatures, breaking CP. This 
condensate is trapped by a barrier in the effective po- 
tential, indicating the presence of metastable states that 
violate CP, which is a relevant point for the possibility of 
observing CP-odd domains in high-energy heavy ion col- 
lisions. We should remark that these metastable states 
were not found in a detailed study of the phase diagram 
defined in terms of the strength of the 't Hooft determi- 
nant, quark masses, temperature and chemical potential 
that has been performed for 9 = 7r within the two-flavor 
NJL model [2B]. The reason for this discrepancy has its 
grounds in the difference in physical content of the two ef- 
fective theories, and is currently under investigation |50j . 

The main effect of a high magnetic background on the 
chiral transition is turning it into a first-order transition, 
even in the case of the CP-odd linear sigma model, deep- 
ening the absolute minimum of the effective potential. 
The behavior of the critical temperature is also affected 
in a non-trivial way, first going up and then dropping for 
larger values of the magnetic field, as was also observed 
in Ref. [SSJ. 

In order to address the physical scenario in the case of 
noncentral heavy ion collisions at RHIC and the LHC, 
one has to incorporate the effects from nonzero chiral 
and baryonic chemical potentials jTU]. It is also crucial 
to describe the nonequilibrium dynamics of the forma- 
tion of the condensates, which will provide the relevant 
time scales and give information on the actual possibility 
of measuring effects coming from the formation of CP- 
odd domains. This issues will be considered in a future 
publication [ST] . 
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